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Abstract 

A recently proposed method of estimating the asymptotic behaviour of QCD 



(N 

D , perturbation theory coefficients is critically reviewed and shown to contain nu- 

'^ ' merous invalid mathematical operations and unsubstantiated assumptions. We 

discuss in detail why this procedure, based solely on renormalization group (RG) 

^ ■ considerations and analyticity constraints, cannot lead to such estimates. We 

j^ . stress the importance of correct renormalization scheme (RS) dependence of any 

meaningful asymptotic estimate and argue that the unambiguous summation of 

QCD perturbation expansions for physical quantities requires information from 

outside of perturbation theory itself. 
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1. Introduction 

Because of computational complications only the lowest two or three terms in 
QCD perturbation expansions are available for phenomenologically interesting 
quantities. Taking into account that the relevant expansion parameter (strong 
coupling constant Us/t^) is of the order of 0.1, and for certain quantities even 
bigger, we face the following two important questions: 

1. how to reliably estimate higher order terms in perturbation expansions of 
physical quantities and, having obtained them, 

2. how to formulate a physically well-motivated summation method in the 
situation when these asymptotic estimates prevent the use of the familiar 
Borel summation technique. 

A number of theorists have struggled over the last two decades to evaluate 
multiloop Feynman diagrams in QCD and other theories and we seem currently 
to be at the border of what can be done with the available analytical as well as 
numerical tools. The state of the art in this field has recently been reached by 
the first three loop QCD calculation of a physical quantity, the familiar i?-ratio 
in e"'"e~ annihilation into hadrons |]I| 

where q^ denotes the square of the center of mass energy in e"'"e~ collisions and Qi 
is the electric charge of the produced quark with flavour i. These, as well as the 
analogous results on certain structure functions sum rules [0 have been used in 
phenomenological analyses of experimental data and the inclusion of NNLO cor- 
rection has been shown to reduce theoretical uncertainties [^ Q and, in the case of 
the Gross-Llewelyn- Smith sum rule, also to improve agreement with experiment 
[0. There seems, however, to be little chance of going in a foreseeable future to 
still higher orders for any of the physically measurable quantities. The quest for 
at least a reliable estimate of these higher order contributions is therefore highly 
commendable, in particular in view of the fact that the relevant perturbation 
series are expected to be factorially divergent with asymptotically sign-definite 
terms. 

In fact the question of the possible divergence of QCD (as well as other field 
theories) perturbation expansions, which goes back to the original argument of 
Dyson [Q, still defies definite and clear answer. The original conjecture that the 
divergence of perturbation expansions is directly related to the discontinuities of 
Green's functions in unphysical domains of the coupling constant has been 
seriously questioned by several authors [0, §, |^. The knowledge of the disconti- 
nuity itself is simply insufficient for the determination of large order behaviour of 




perturbation theory coefficients and one needs much more detailed information 
on the behaviour of Green's functions in the vicinity of the corresponding cut in 
order to estabhsh this relation. 

We share the viewpoint of Stevenson ^ that there is little evidence that QCD 
perturbation expansions for physical quantities do indeed diverge. The examples 
discussed in |^ indicate that in fact most of the information on the analytic 
properties of relevant Green's functions near g'^ = 0, |lT| may actually be 



invisible in perturbation expansions themselves. Moreover, the fact, contained 
already in and later rediscovered and advocated by Stevenson [§], that the 
summability of perturbation expansions is not the same as the possibility of 
recovering from such a series the corresponding "full" physical quantity, has now 
gained wider acceptance. As the summability of QCD perturbation expansions 
is determined by the large order behaviour of their coefficients we primarily need 
to know which kind of information is necessary and sufficient to derive it. 

Quite recently two new attempts have been undertaken in order to derive the 
asymptotic behaviour of the perturbation-theory coefficients for certain physical 
quantities p, [1^. They use completely different ideas but both arrive at the 
same conclusion: QCD perturbation expansions are indeed factorially divergent. 
The novel argument of is tailored for the special case of particle fractions in 
e"'"e~ annihilation and can not be used for such simple quantity as (|l]). Contrary 
to the original approach of it, however, doesn't employ any information on 
analytic properties of corresponding Green's functions. Although it is based on 
several strong and questionable assumptions, which render its conclusions open 
to doubt, we can at least imagine that more sophisticated procedure can be 
formulated along similar lines. 

This is not the case with the procedure suggested in ||12|, which is claimed to 
answer both of the above questions. In a subsequent paper JTsf, West uses his 



estimates to formulate a summation procedure with truly remarkable properties. 
In particular they make the use of a few lowest terms in (|l]) a reliable approxi- 
mation even in the low q^ region where the expansion parameter as(g^)/7r is so 
large that this approximation is normally considered meaningless. 

Soon after it appeared, |]T2| was criticized by Brown and Yaffe |T^, who have 



argued that the assumptions employed by West, i.e. the analyticity constraints 
coupled with RG considerations, allow one to find the relation between the Borel 
transform of (0) and that of the dispersive part of n(g^) but are not sufficient to 
unambiguously determine the asymptotic behaviour of the coefficients r^ them- 
selves. We fully agree with these authors but consider it useful and illuminating 
to go through the procedure of |jl2| carefully step by step in order to see exactly 
where and why it breaks down if carried out correctly. This will be done in Sec- 
tion 3. Our analysis will also show why there is no hope to remedy the derivation 
in |T2|. Finally in Section 4 we shall comment on the summation procedure sug- 
gested in [|13| , which uses the factorially divergent nature of the coefficents r^ but 



is in principle indepedent of the particular way of deriving it. We shall argue 



that this algorithm is just one of an infinite number of equally plausible "sums" 
of (|l|) compatible with the factorial growth of r^. 

In renormalized field theories like QCD the question of the large order be- 
haviour of the coefficients like r^ in (|lD is complicated by their inherent RS 
dependence. Consequently even the convergence property of the expansion (|1|) 
does in fact crucially depend on the RS chosen. We therefore start our analysis 
with general remarks concerning the constraints following from the formal RG 
invariance of expansions like (|l]) which any asymptotic estimate of r^ must obey 
in order to make sense. 

2. RG constraints and RS dependence of asymptotic es- 
timates 

The expansion parameter a^ as well as the coefficients r^ in (0) depend on the 
RS chosen. A part of this dependence is labelled by the scale fj, but the full RS 
ambiguity is much wider. As this point is often a source of misunderstanding let 
us specify very carefully what is understood under this term and where it differs 
from the so called renormalization convention (RC), the notion introduced in 



1^ . We have decided to include this detailed exposition after getting acquainted 
with West's response [^ to |T^ and other critics |1^ . 
We start with the usual RG equation 

'^^^^'^ - m = -g' (bi + hg' + hg' + ■■■), (2) 



dln/i 
which represents basically a definition of the expansion parameter 

a{^?) = a.(/i2)/7r = g\fi')/An^ 

used in (|l]). While 6i as well as 62 are unique numbers, fixed by the number Uf of 
quark fiavours (we consider massless quarks only), all the higher order coefficients 
bn, n > 2, are arbitrary numbers defining the renormalization convention. For 
physical quantities the notion of the renormalization scheme involves, besides the 
arbitrariness in the /3-function coefficients, another degree of freedom, not related 
to them. To exemplify this point let us just recall the difference between the MS 
and MS RS. Both RS have exactly the same /3-function coefficients but they 
are nevertheless associated with different values of the couplant a^fi"^) as well as 
the coefficients rk,k > 2. In fact there is an infinite set of RS which have the 
same /3-function as MS, but which correspond to different values of rfc,fc > 2. 
This dependence of the coefficients r^ on the RS chosen is formally to the k-th 
order compensated by that of the couplant a(/i^, RS) in a way which guarantees 
the internal consistency of perturbation theory |T^ . There are various ways of 
parametrizing this degree of freedom, the most conventional one using just the 
value of renormalization scale /i. 



Mathematically, this additional degree of freedom in the definition of the RS 
is related to the simple fact that the RG equation (^ defining the renormalized 
coupling constant 5'(/i^) (and thus ^^(yU.^)) has, even for fixed values of all the 
coefficients 6j, an infinite number of solutions. For any well-defined r.h.s. of @, 
like for instance in the so called 't Hooft's RC, (or in any of the "analytical" RC 



discussed in [0) where all free higher order coefficients &«, i > 3 are set identical 
to zero |T^, the behaviour of g{fi^) for fi ^ oo is essentially fixed by the first 
term in (Q). At finite /i, however, the differences between various solutions of (|^) 
may be arbitrarily large. 

Specifying the RS then means besides fixing the /^-function coefficients also 
choosing the renormalization scale /x. Variation of as(/i^) with this scale is then 
compensated to any finite order by that of the coefficients r^. For a fixed /i there 
is, however, an infinite set of the values g{n'^) corresponding to different solutions 
of (^. So even setting fi equal to some external kinematical variable (like the 
total CMS energy y/q^ in ([^)), we still face the necessity to specify which of the 
solutions to (^ we have in mind. This fact is frequently overlooked and one often 
finds statements that the fixing of the scale n and the /3-function coefficients b^ 
defines uniquely the renormalization scheme. This is, however, inaccurate and the 
selection of the "referential" RS, i.e. the selection of one of the solutions to (0) as 
referential is as important as the choice of the scale /z itself. Indeed, the variation 
of the scale /i for fixed referential RS (RRS) is equivalent to the variation of the 



RRS for fixed /i. In the approach suggested and advocated in ||T5[ , we select once 
for all one of these solutions and label different RS by different values of /i and bj. 
The selection of this referential RS is of course merely a matter of bookkeeping 
and has nothing to do with the actual choice of a particular RS. The RS is thus 
uniquely defined by either of the sets {/i, bi, RRS}, {a, 6j} or {r2, bi}. In the last 
two cases there is no need to select the RRS. 

There are various ways how to label the solutions of (0), the most convenient 
one employing the renormalization group invariant, dimensional parameter de- 
noted usually A. Each of the solutions of (^ is then associated with one particular 
value of A (we denote it Arrs) and the couplant a{fi'^) is then in fact a function of 
the ratio /i^/A|^j^g. The precise definition of Arrs being again a matter of conven- 
tion, we can define it, for example, by the relation al^"^ = A^Rg) = oo. Knowing 
the dependence of the couplant a(/i^/A^; RRS) on both the scale and the RRS 
we come now to that of the coefficients r^. For them the internal consistency of 



perturbation theory implies the following recurrence relations [15 
drfc(gV/i2;RRS) 



dln/i 



8TT%{k - l)rfc_i + 32TT%2{k - 2)rk-2 + 0(63)- (3) 



For r2 this yields (supplemented with the initial condition tq = ri = 1) the 
following explicit dependence on /i 

r2(gV/^'; RRS) = 47r26iln(/iVg2)+^2(l;RRS) (4) 



= 87r%ln(^/ARRs)-pi(g^), 

while for r^ we get 

r3 = p2-C2 + (rs + c/2)2, (5) 

where c = &2/^i and C2 = 63/&1. As a consequence of (^ analogous algebraic 
relations hold also for higher order coefficients rfc. Similar consistency conditions 
can also be obtained for the variation of r^ with respect to other RS parameters, 
i.e., bi,i > 3. 

The quantities pi and p2 (and analogous quantities at higher orders) entering 
the formulae (^ and (^ are RS invariants, pi being a function of q^ while p2 (and 
all higher order invariants pi,i > 2 as well) is a pure number. As fi and r2 are 
(in a given RRS) in one-to-one correspondence we can, instead of /i, alternatively 
use the value of r2 itself to label the same degree of freedom as fi. In other words, 
the fixing of r2 is an alternative way of fixing (together with the coefficients 6„, 
n > 2) the RS. We stress that it is already the value of r2 which is different in 
the MS and the MS RS, although the corresponding /^-functions have exactly the 
same coefficients bi. 

We mention this alternative labelling of the RS by means of {r2,bi} for the 
following reason. It is in fact difficult to translate the results on the coefficients r^ 
obtained by means of any method not explicitly using Feynman diagram renor- 
malization (like the one in [Q) into the language of the conventional perturbation 
theory and vice versa. For instance, if, in West's approach, one were to carry 
out the evaluation of r^ in, say, MS RS there seems to be no obvious meaning of 
the term "minimum subtraction" , which is inherent to the usual Feynman rules 
calculations. The same holds of course for any other conventional RS. Using (§) 
we can, in any selected RRS, define /i associated with the MS RS by demanding 
that the corresponding r2(/i^/g^; RRS) = 1.41 but this is possible only if the de- 
pendence of r2(/i^/g^; RRS) on /i is consistent with (|^). In the notation where fj, 
is fixed, being set equal to, say, y/q^, similar statements hold for the dependence 
on the RRS. In [|I^, however, the asymptotic behaviour of the coefficients r^ 



TT ^ ' in + b' 



rn{l) - --^{^r.'eb,r~' ;A!iZ^, (6) 



is manifestly RRS independent. Indeed it was the striking lack of any RRS 
dependence of (P) which has originally attracted our attention and cast doubts 
on their validity. The absence of appropriate RRS dependence of r„(l), prevents, 
however, the association of (^ with any well-defined RS. This crucial fact has 
not been properly appreciated in [|T^ . 

To see clearly the implications of the consistency conditions (|^) for the RRS 
dependence of the large order behaviour of the coefficients r^ let us consider the 
(for technical reasons somewhat simplified) case when 6^ = for all A; > 1, i.e.. 



also 62 = 0. Then (^ can be explicitly integrated to yield 

j=o \ 3 J 

where r = 87r^6i ln(/i//io), /io is some "initial" value of /i (corresponding to some 
RRSo) and rk{q'^/f^l) are the associated expansion coefficients of (^ in RRSq. 
The above formula can simply be rewritten to express the RRS dependence for 
fixed ix^ = q^: 

rfc(l; RRS) = ^ r^rfc_,(l; RRSq) (^~^\ (8) 

i=o \ 3 J 

where now 

r = 87r26iln(ARRSo/ARRs). (9) 

Assuming, as an example, the factorial behaviour of the coefficients, rfc(l; RRSq) 
fc!, we get from (^) for general RRS 

rfcfl; RRS) = k\y ^ => lim ^ = expr. (10) 

^ %] k^oo fc! 

This means that although the character of the large order behaviour of the coef- 
ficients rfc(l;RRS) is not changed (in the above example it is merely multiplied 
by the fc-independent factor exp r) when we vary the RRS, the magnitude of the 
coefficients rfc(l; RRS) obviously is. As r can in principle be any real number, the 
magnitude of the coefficients rfc(l; RRS) can, in our example, be arbitrarily large. 
Although the RRS dependence of the large order behaviour of ^^(l) is irrelevant 
for considerations concerning the summability of our series, it is crucial when 
we are interested in their numerical values. For the realistic case with &2 7^ 
the dependence of rfc(l; RRS) on the RRS is more complicated than ( p!0D but its 
main feature i.e. the arbitrariness of the magnitude of rfc(l;RRS) on the RRS 
persists. We want to emphasize this point very strongly as in his reply [^ West 
explicitly claims that the large order behaviour of the coefficients ^^(1) should be 
"scheme invariant" with "scheme dependence" entering only via the dependence 
of nonleading corrections on the free coefficients bk,k > 3. The relation (|10|) 
demonstrates that this claim is definitely wrong as the dependence of rfc(l; RRS) 
on the RRS may change the magnitude of them essentially at will. For example, 
the value of r corresponding to the change from the MS to the MS RRS equals, 
for Uf = 5 

T = %Ti%i ln(Aj^/AMs) = 3.83 In 2.65 = 3.73. (11) 

There is thus no obvious reason to compare r^ = —13.4 from @ with its exact 
value in the MS RS (where r3(l; MS) = —12.8) and not, for instance, in MS which 
has exactly the same /3-function coefficients but where, however, r3(l;MS) = 
16.3. The agreement of (H) with r3(l;MS) should therefore be considered as a 



pure coincidence. This has also been demonstrated in |T^ on the basis of the 
comparison of the corresponding Uf dependences of (H) and the exact resuhs 0. 
There is in fact no agreement of (D even in the MS RS and even following 
uncritically the whole derivation of |T^. The "correct" (i.e. disregarding all the 
objections we shall discuss in the next section) result thereof reads 

rM - -^(4-=^M"-' ^^, (12) 

which differs from (^ asymptotically by a factor {n^^ )/2 and by a factor of 
about 1/4 for n = 3. So even within the framework of West's approach the 
starting point for his considerations, i.e., the coincidence with exact results of 
is actually absent. On the other hand it is clear that for any value fs there 
is always such a RRS that r3(l;RRS) = r^. We also stress that the ambiguity 
in association of (^ with a particular RS has nothing to do with the neglected 
higher order /3-function coefficients bk but reflects the fact that even for a given 
RC, i.e. for fixed fi and bk,k > 3 there is still a degree of freedom connected with 
RRS. 

In perturbation theory the quantity (^ is described by a pair of expansions 
(|l]) and (H), neither of which has a physical meaning of its own. Except for 
the first two coefficients in both (|I|) and (0) all other ones are arbitrary, subject 
only to the internal consistency constraints like (|^) which follow from (|^). We 
can choose the RS={/i, 6j,RRS} in such a way that either (|I|) or (|^) is a well- 
defined, convergent series but we can not do this for both of them simultaneously. 
For instance, the "effective charges" approach of Grunberg [|l^] sets r^ = for all 
A; > 1 so that one has to worry about the convergence property of (^ only. On the 
contrary, in 't Hooft's RC 6fc = 0, k > 2 so that the definition equation (j^) of the 
couplant a{fi'^) is well-defined and the nontrivial question concerns the expansion 
coefficients r^. In a computationally simple RS, like MS, both of the expansions 
(|I]) and (0) are presumably divergent. This makes their use for the purpose 
of summation considerations rather unsuitable as the couplant a^/vr itself is ill- 
defined in the limit of infinite order. Consequently, any statement concerning 
the large order behaviour of the coefficients r^ is inextricably connected with 
that of the coefficients bk of (|^) and makes sense only when it is clearly and 
unambiguously related to some definite RS. 

3. Discussion of West's approach 

We now turn to a discussion of the West approach advocated in [l^. (For reader's 
convenience we basically use the original notation of fT^.) We start from the 
definition of the vacuum polarization function Il{q'^,g'^) in terms of the time- 
ordered product of currents: 

{q^.q. - 9^..q') ^(g^/) = ijdx e"^^ (0|Tj^(x)j,(0)|0), (13) 



where j^(x) is the electromagnetic current of the quarks in massless QCD. The 
function 11 is assumed to be analytic in the whole q^ complex plane cut along the 
positive real axis. The imaginary part of 11 is related to the ratio R (see eq. (|l])) 

Im n(gV/i', as(/i')) = -^R{qV ^\ as(/x')). (14) 

We further introduce the function 

where t = ln(g^//i^). R and D satisfy the homogeneous renormalization group 
equation 

[fi^ + P{g)^]R{qVl^',as) = (15) 

and analogously for D, where the Callan-Symanzik function (3 is defined in (0). 
As a consequence the functions R and D depend (in massless QCD) on the single 
variable z = (g^//i^)e^-^*^^^ |]TD|, |T2|. Here the function K{g) is defined by 

r9 rln' 1 h' 

K{g)= 7^^;^r^-^M&i/^' + &2)+---, (16) 

Jga l3{g') 2hig^ 2 

where h' = 62/^1- (In obtaining eq. (|l^) a convenient choice oi go has been made.) 
It is useful to define the running coupling constant g(t) by 2K{g{t)) —2K{g) = t. 
Then the eq. (|l|) can be rewritten as 

00 
i?(gV/x^a,(^2)) = Rihasiq')) ^ (3Ea') E ^-(1) («.(?')A)"- (17) 

i n=0 

The dispersion relation for D can be given the form 



fOC 

D{qyf^\g') = A{qyf,\g') dz 

JO \z — 



A{q^lli\g^T " ' 

where 

J{z) = -^^R{q^l^\a,{^i^)\ (19) 

and 

A{q^l^\g^)=^-^^^^^\ (20) 

(In the following we shall use the reduced notation A(A;) = A[cj^ j ^ ^ 1/^); where 
k = l/g\) 

The aim of West's considerations is to use eq. (|T^) in the analysis of large n 
behaviour of the expansion coefficients rf„ defined by formal (asymptotic) expan- 
sion 

00 

D{qVf^\ /) ^ E (-l)"^n(gV/i') /"• (21) 

n=0 



The expansion coefficients r„, defined in (|T]), are then connected with the dn by 
the relation 



(E^n r„ = -(-47r2r+i^[Imrf„+i + ^Imd„ + ---], (22) 

\i J T^nhi hi 

as a consequence of the equation 

This suggests that we can get the large order behaviour of r„ once we know that 
of dn (see however |]13[)- 

We now focus on three crucial points in the method of calculation of the 
coefficients dn- 

(A) If eq. (pTD can be understood as an asymptotic expansion of D then the 
dn are given by (see Wightman |1T0| ) 



dn = - lim f Jf:i-g^ + g^,r-Diqy^^^g^ (24) 

30^0+ JC(go) Zm 

C{go) is a closed contour around the point gQ. The size of the contour has to be 
chosen sufficiently small in order not to hit the complicated singularity structure 
of the function D around the origin. In particular, C{go) must not leave, with 
go approaching zero, the wedge bisected by the real axis and bounded above and 



below by circles that are tangent to the real axis at the origin|jTO|. (There are 
arguments that this horn-shaped region is the analyticity region of the function 
D [0.) Making the transformation k = l/g"^ and choosing the corresponding 
contour C'iko) around the point k^ = Ijg^ (see Fig.l), we can rewrite eq. (^) 
into the form 

dn = lim / 1^ r-1 A;o"+^ {k - fco)-^-" D{q^ / ^^\ l/k). (25) 

fco^oo Jc'(fco) 27r« 

One particularly convenient choice of C'{ko) is the contour Cq (see Fig.l) which 
is independent of the position of ko provided that ko is sufficiently large. (We 
denote as Cq the corresponding contour in the ^f^-plane.) We do not extend the 
contour Cq up to Re A; — » — oo, in order to avoid the problem with a singularity of 
D a.t \k\ = oo. Thus, in the contour C[ depicted in Fig.l one cannot neglect the 
contribution along the direction perpendicular to the real axis at Re /c —^ — cxo. 



The contour mentioned in |16] should therefore be understood as the Re /c — *■ — oo 
limit of C(. 

Obviously in eqs.(p4D and (|25|) one can not interchange the order of the limit 
kf) —>■ oo and the integration. Otherwise we would obtain the following result 

dn=<f ^i-kr-'Dig'/ftM/k). (26) 







10 



Using the expansion ( pT] ) in the neighbourhood of Re /c ^ — oo we conclude that 



the integral (^) does not exist for any n > 0. Rewriting (p6D back to the g^-plane 
we obtain the expression 

dn = ~f ^(-/)— ^D(gV/^^/), (27) 

JCo ^TCl 

which coincides with West's formula (with the integration contour properly de- 
fined). Thus, the dn presented in |jl2| cannot be understood as well-defined 
quantities. 

(B) The next step in West's approach is the substitution of (|18|) into ( pS]) 



and the interchange of the order of integrations. We shall however show that the 
dn computed in this way are incorrect (this suspicion was also voiced in ref . Jl^ ) • 
For this purpose we introduce the quantity 



d 
dn = lim / dz t [z) 

ko^oo Jo 

where 



POO 

dn= lim / dz f{z) — ^{n,t,z,ko), (29) 

fco^oo Jo at 

$(n, t, z, ko) = <p dk -^— -— . (30) 

^ ' ' ' ^ 27ri Jc' (k - fcn "+^ z - A(k) ^ ' 



.„ {k - koY+^ z - A{k) 
Note that the formal limit /cq — > oo would recast this relation to the form 



*("•'--' -2^ 4*1^' p^' 



used in pf. 

Applying the Cauchy theorem to (pOf) we get 



C 1 fin i-"-! 
<l>(n,t,^,M = k-+'\- ^ 



n\ 9A;" z-A{k) 



k=kQ 



k;-\hkp + b2) 



Sc (32) 



{kp - kor+^Aik,) 

= ^i{n,t,z,ko) + ^2{n,t,z,ko). 

The contour is chosen such that at most one of the poles given by 

A{k) = z- e*+2^('=) = 0, (33) 

can lie inside it. We will denote its position by kp (see Appendix). 6c' is a 
contour dependent factor, 6c' = 1 if the pole is inside C, otherwise 6c' = 0. 

11 



This contour- dependence is apparently undesirable and it is the first signal that 
dn given by (|29|) differs from dn {dn is contour-independent provided that ko lies 



inside Cq) 



Using (|3^) it is easy to show by induction over n that 

d 



lim / dz f{z)—^i{n,t,z,ko) = dn. (34) 



Thus, we get 



dn - dn= lim / dz f{z) — $2(n, t, z, ko). (35) 

fco^oo Jo at 

We conclude that this contour-dependent difference is given by the pole contri- 
bution to dn- The position of the pole depends on the value of z. (When 62 ~ 
we have kp ^ bi{lnz — t).) In general, the pole moves to the right in the k- 
plane when z increases. The integration in ( |53| ) can therefore be restricted to an 
interval [zi, Z2) defined by the occurrence of pole inside C . In the case of the 
contour Cq depicted in Fig.l we have Z2 = C)0. For our aims it is sufficient to 
prove that the difference (^) does not vanish. In fact we prove in the appendix 
that the difference can be divergent for some contours. Consequently, one cannot 
interchange the order of integrations in (pS]). The same objection has been raised 



in |T^ (in this paper, however, the contribution of the pole to dn is ignored). 
(C) In |T^ the integral (^) has been estimated by the saddle point technique 



which might be expected to provide a large-n approximation to (|3^). An essential 
feature within the saddle point method is a convenient choice of the integration 
contour. Representing the integrand in the form of exph{k), where h{k) is a 
complex function analytic in some vicinity of the contour, we should choose the 
contour so that the following conditions are satisfied: 

(i) The integration contour can be deformed (without changing the value of the 
integral) so that it passes through the saddle point kg. 

(ii) Re h{k) has its maximum on the contour just at the saddle point kg. 

(iii) The contour can not be deformed so that a new maximum of Re h{k) is 
below Re h{ks). 



A detailed inspection of the integrand in (31) reveals the fact that the con- 
dition (ii) is not satisfied if we integrate in the direction perpendicular to the 
real axis in the complex fc-plane. The curve of the steepest descent passes along 



the real axis. This fact invalidates the approach of [|T2| where the integration 



contour is not specified and seems to be chosen perpendicular to the real axis at 
the saddle point. Moreover, a detailed analysis of the saddle points shows their 
strong dependence on the value of z. This crucial fact was also overlooked in [0, 
where k^ ^ hi{n — 1 +6') was taken. We illustrate the situation in Figs. 2-5, where 
the positions of poles and saddle points of the integrand in ( piD are displayed for 
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different values of z and for n = 4. One can notice an interesting interplay of 
poles and saddle points resulting in the situation when two complex conjugate 
saddle points become leading {i.e. with maximal Re /c) for sufficiently large z. 
The z dependence of the saddle point kg = ks{z) does not allow one to choose 
any universal (^-independent) integration contour in the /c-plane. Hence, the use 
of the saddle point method is plagued with the occurrence of the moving leading 
saddle point ks{z). Consequently, one cannot obtain the relation ($(fc) is defined 



inlll) 

1 r 2 1'/' 

dniqVl^')^- ^^^^ k:-'Diqyf^\l/k,)cosnn, (36) 

which plays a crucial role in the derivation of West's result. (We have added 
the factor 1/2 which is missing in |l^.) It was just this relation which made 
it possible, employing (|2^ ) and (0), to derive the simple relation (|) for the 
coefficients of the expansion (|l^. We would like to stress that even within West's 
approach, based on the saddle point method, the (i„ will also get contributions 
that cannot simply be related to the value of D(g^//i^, 1/k) for some large k (in 
( P^D kg ^ oo when n — »• oo). This follows from the 2;-dependence of the leading 
saddle point kg mentioned above. The same argument can also be used for the 
case of finite ko in (|30|) . 

Moreover, the pole of the integrand in (^) generates, at finite /c = /cq, an 
additional saddle point which disappears when fco — >■ oo, and for some z, becomes 
the leading one. For these values of z saddle point method can be used to estimate 
the integral (|30|) because the conditions (i)-(ii) are satisfied. This is another signal 
that the saddle point method of approximating the dn is not applicable in the 
form advocated by West. 

4. Summing perturbation series in QCD 

Using the asymptotic estimate (H) West has in his subsequent paper suggested 



rH| a way of handling the divergence of QCD perturbation series which has two 



remarkable properties: 

1. it is extraordinarily close to the "correct" sum using only a few lowest order 
terms and, moreover, 

2. the number of terms to be taken into account in his procedure is inversely 
proportional to the value of the couplant, which implies that it decreases 
as we approach the infrared region, i.e. as g^ — > in (2). 

These properties are remarkable and, were they true, would have serious reper- 
cussions as the conventional wisdom tells us that at low g^ the contribution of 
nonperturbative effects are expected to become dominant [jl9|, ^ . As the message 
of |TB[ is, from the viewpoint of phenomenological applications, very appealing we 
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wish to point out several simple and obvious facts which show the arbitrariness of 
West's construction. Basically what he does in |T^ is that he inadmissibly gener- 



alizes the conclusions valid for one particular definition of the sum, which really 
has the above-mentioned properties, to the case of QCD perturbation expansion 
(|1|). Let us recall the basic fact that for the divergent asymptotic series 

oo 

E «n-^" (37) 

n=0 

there is no "natural", unique, sum. There is an infinite number of functions F{z) 
which have the series (|37D as its asymptotic one and these functions while sharing 
the same behaviour in the vicinity oi z = may for any fixed z 7^ be arbitrarily 
far apart. Perturbation theory itself does not tell us which one of them to choose 



and we have to look somewhere else for ideas what to do with it pp] , |2T| . Without 
first defining what we mean under the "sum" of the asymptotic series (^), there 
is no sense in talking about the "error" of any finite order approximation. To 



claim, as is done in 131, that this error, defined as 



^■5 



R 



N 



iV- 



z - 



E ««^" 



n=0 



(3^ 



can be estimated by 



i?^ ^ a^z^ (39) 



is therefore unsubstantiated. 

The asymptoticity of the series (|37| ) merely means that 



R^(z) 
lim ^ ^J_, = (40) 

but, for a fixed z ^ 0, R^ may be quite arbitrary. In other words any statement 
concerning the error estimate (|38|) is meaningless unless we unambiguously specify 
the definition of the sum F{z) appearing in ( ^81) as well. This holds equally well 
for Borel summable as well as non Borel summable series. In [^ the author 
bases his claims essentially on the property of a particular Borel summable series 
with the coefficients r^ = {—l)'^k\. He demonstrates that by choosing the number 
N{a) of terms in (|38D , where this value is roughly inversely proportional to the 
expansion parameter a and given by the local minimum of the contribution of anz"' 
to the partial sum, he gets a very good approximation of the corresponding Borel 
sum. This is by no means the only possible approach, even for Borel summable 
series: one needs further information on the function F{z) itself (i.e. information 
not contained in the perturbation expansion) to guarantee the uniqueness of the 
Borel sum. To generalize the properties of this toy example series (coupled with 
one ad hoc definition of the sum) to the non Borel summable series in realistic 
QCD is an entirely arbitrary procedure. Moreover, the procedure suggested in 
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IH] has numerous concrete shortcomings. We mention just the most important 



ones: 

• It is obvious that by adding to this toy example series any convergent one 
we get the expansion with the coefficients which have the same large-n be- 
haviour while to any finite order the partial sums may be arbitrarily far 
apart. Any meaningful statement concerning the accuracy of a given sum- 
mation technique must take into account the true values of the expansion 
coefficients, not merely their asymptotic behaviour. To determine the value 
of N{a) from this asymptotic behaviour is unsubstantiated. 

• The simple picture in which the contributions of the terms ttnZ^ first mono- 
tonously decrease and then starting at N{a) increase forever is a property 
of the toy example mentioned but in realistic QCD the behaviour of this 
contribution may be much more complicated. In fact, for a function F{z) 
to have ( |37D as its asymptotic series we can take N{a) entirely arbitrary, 
subject only to the condition that as a — ;> 0, N{a) — >• cx). For a ^ 
all these procedures do, of course, converge to the same result but for any 
nonvanishing a they may yield vastly different numbers. Consequently, also 
the estimate (17) in [T^ of the error of the procedure suggested is valid for 
the example discussed but has no general validity. To give his claims some 
justification West would have to argue why his definition of the estimate is 
to be preferred to the infinite number of other possibilities. 

• For the asymptotic series of the type r^ ~ /c!, expected to occur in QCD, 
there is no "natural" sum like the Borel one for the sign alternating series 
Tfc K, [—l)^k\. To define in these circumstances the estimate by means of the 
same, above mentioned, procedure is unsubstantiated and has no physical 
motivation behind it. 

• Even if all the coefficients r^ in (2) were known exactly (as in the toy exam- 
ple of the series with coefficients r^ = (— l)'^fc!), the summation procedure 
suggested in JlSl would still be itself badly ambiguous as the value of N{a) 
and consequently also of the estimate defined as 

N{a) 

E ru{as/7^f (41) 

fc=0 

would depend on the RS in that the truncation procedure is carried out. 

5. Summary and conclusion 

We have examined in detail the method recently proposed by West for the de- 
termination of the large order behaviour of perturbation expansion coefficients in 
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QCD. Going carefully through his derivation, which is based exclusively on the 
renormalization group considerations and analyticity properties of the vacuum 
polarization function, we have identified several points where strong and largely 
unsubstantiated assumptions have to be made in order to arrive at the result (P); 
in particular 

(A) the interchange of the limit fco — ^ oo and the integration over k in (p5|); 

(B) the interchange of integrations in the double integral in (^), which results 

in the inadmissible contour dependence of the result; 

(C) the use of the saddle-point technique for the estimate of the contour integrals, 

in spite of the ignorance of the 2;-dependence of the positions of the leading 
saddle points. 

A first signal that there is something wrong with the result (|^) comes already 
from the fact that this formula for r„ is independent of the RRS chosen. However, 
as is obvious from his reply |]l6l to the early critical voices 0, , West considers 
this property of his formula (^) on the contrary a merit. To clarify this crucial 
point we have discussed at length in particular the correct RS dependence of 
any consistent large order behaviour of the coefficients r^. In general, we thus 



agree with the claim of Brown and Yaffe[0] that the information exploited by 
West is in principle insufficient for the derivation of the asymptotic estimates of 
these coefficients. Compared with their work, which shows what can be derived 
from these assumptions, we have attempted to demonstrate why the construction 
suggested by West doesn't in fact work if carried out with proper care, and why 
there is no conceivable way to improve it and still arrive at an asymptotic estimate 
of the type (^. The large order behaviour of QCD perturbation expansions 
thus remains to large extent an open question which certainly deserves further 
attention. In our view, to give these expansions good mathematical as well as 
physical meaning we have to go beyond the purely perturbative framework and 
look, as in |2y, ^, ^ , to their interplay with various nonperturbative effects. 
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Appendix 

First, let us derive the equation for the boundary of the principal strip in the 
/c-plane. The relation 
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transforms, at g^ fixed, the complex z-plane into a strip in the fc-plane. K{g) is 
given by (|16]). Let us introduce the notation 

k = u + iv, g7/i2_ae^'^, lna = t, z = re'^ (43) 

and take the logarithm of ( ^2] ) (In stands for the principal branch): 

In a + u/bi - -b' ln[{biu + 62)^ + blv^] = In r (44) 

and 

I-kN + v3 - ^ + w/61 - b'arg{bik + 62) = (45) 

The borders of the strip are given by the inequality 

- TT < ^ < vr, (46) 

which represents the whole complex cut z-plane. Choosing A^ = and </? = 0, we 
obtain the implicit equation for the upper bound of the strip in the form 

V — — arctan(- — ) — nbi -n6{~biu — 62) = 0, v > 0, (47) 

Oi biu + 62 bi 

the lower border being symmetric with respect to the real axis. 

We now turn to the contour dependence of the difference defined in (p5[) : 



rz2 Q 

Adn = dn-dn= lim / dz f{z) — ^2{n,t,z,ko), (48) 

ko^ooJzi at 

where 

and for z G (-21,-22) the pole kp occurs inside the integration contour. We note 
that a usual regularization is understood in the denominator factor [kp — k^) of 
(^91) : we can add a small imaginary part to k^. The function kp{z) solves the 
equation (|33D. Denoting kp = Up + ivp one can show that the positions of poles 
satisfy the equation 

-\\nvl + \ \n[sin\vp + 7)] - 5 = 0, (50) 



tan({;p + 7) 2^2 



where 



Vp = 


hi 




7 = 


= 27cNbl/b2, N = 


= 0, 


B = 


52 

= l + \nb2 + T^{\nz- 
02 


-t) 
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provided that f p + 7 7^ mr, n = 0, ±1, ±2 . . .. Knowing Vp the real part Up can 
then be calculated from the equation 

^^""^| = tan(';+7)-' ^''^ 

A different regime can be observed if Wp + 7 = nir. These are the relevant poles 
because the integration contour is situated in the right half plane. (We consider 
only such contours which can be hit only by one pole, i.e., which are placed in 
the strip.) The poles on the real axis satisfy the equation 

Up-^ ln[(l + UpY] -B = Q (52) 



Here the corresponding Vp equals, of course, zero. In the 5-flavour case, i.e., 
61 = 0.0485 and 62 = 0.00155, we get that for z > 71.024 the equation (^ has 
two solutions, one with Up positive, the other one with Up negative. They move 
to opposite sides if z increases starting from the origin. If 2; < 71.024 we get no 
real poles in the strip. However, one can find two complex conjugate poles in 
the left half of the complex /c-plane. The ^-evolution of the poles is depicted in 
Figs. 2-5. As is seen there are two poles which collide ai z = 71.024 from either 
side along the imaginary axis and scatter along the real axis. In the following 
we will consider the pole moving along the positive real axis; let us denote its 
position by kp{z), for z > 71.024. 

Taking the derivative of (^Q]), substituting the variable y = kp{z) into the 
integral over z and defining 



2/2 _ y"- 



Fn,m{h^yi,y2) = dy f{y)- — - (53) 



we get 



Adn = , lim k^+^ {bi{n + l)F„_i,„+2 + b^n + l)F„_2,„+2 

- bl{n - l)i^„_2,n+l - &i&2(n - 2)F„_3,„+l}. (54) 

Here yi = kp{zi), i = 1,2 are the intersections of the path of the pole with the 
contour, f{y) = f{kp^{y)); for simphcity, we have omitted the arguments of the 
functions F„ „j- 

Obviously, the behaviour of Adn depends on the contour chosen. The general 
proof of the behaviour of Adn is, however, complicated because the function f{z) 
is unknown. In realistic situations f{z) is positive and bounded from above by a 
constant. In the asymptotic region we consider the usual asymptotic expansion 
([TTD; which suggests to expand f{y) in powers of 1/y. It is therefore instructive 



to consider the case f{z) = l/yK For / = 0, we have 



F \- - - V 



"0 



^ ^ "o vjy^-j 



X 



^ [{y2 - ko)""-^-^ {yi - ko)^-^-^ J ' 

We can now check the different contours; simple examples are (a) yi — ko = const 
and y2 — kf) = const and (b) the contour Cq, i.e., yi = const and 1/2 = cxd. It 
can be seen from eq. (|^) that in these examples Adn diverges when k^ —>■ 00. 
A similar analysis can also be carried out for f{z) = l/y\ I > 0. The divergent 
parts cannot compensate each other, so we are led to the conclusion that the Adn 
are typically divergent. Concluding this appendix we would like to stress that 
the contour dependence of dn is quite unacceptable. 
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Figure Captions 

Fig.l The integration contours in the /c-plane (see Sec. 3). 



Fig. 2 The 2;-evolution of the positions of poles inside the principal strip in the k- 
plane. Two complex conjugate poles "collide" at the origin (for z = 71.024) 
and "scatter" along the real axis. 



Fig. 3 The 2;-evolution of the positions of the leading saddle points. Two real 
saddle points "collide" and set out to the complex plane. The saddle point 
near k^ completely disappears for k^ -^ oo (see eq. (0)). The stars corre- 
spond to the positions of the saddle points at z = 71.024. 



Fig. 4 The plot of the real part of the logarithm of the integrand in eq. (pOl). 
The values chosen are: n = 4, Inz = 10, t = and k^ = 1. 



Fig. 5 The contour plot corresponding to the situation of Fig. 4. The resolution 
is not sufficient to detect the pole located to the left of the origin. 
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